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Abstract 

The production of jets, and high-pr leading pions from jets, in d-Au collisions at the Relativistic Heavy-Ion Collider 
(RHIC) and p-Pb collisions at the Large Hadron Collider (LHC) are studied. Using a modified version of the event 
generator PYTHIA, in conjunction with a nuclear Glauber Monte-Carlo event generator, we demonstrate how events 
with a hard jet may be simulated, in such a way that the parton distribution function of the projectile nucleon is frozen 
during its interaction with the extended nucleus. Using our hybrid Monte-Carlo event generator, we demonstrate 
that the enhancement in RpA seen in peripheral events at RHIC and at LHC, as well as the depletion in central or 
semi-central events, is mainly due to “mis-binning” of central and semi-central events with a jet, as peripheral events. 
This occurs due to the reduction of soft particle production caused by a depletion of energy available in a nucleon 
(of the deuteron in the case of d-Au collisions), after the production of a hard jet. This represents a form of “color 
transparency” of the projectile nucleon, which has fluctuated to a state with fewer and harder partons, in events which 
lead to jet production. We conclude with discussions of the form of multi-parton correlations in a nucleon which may 
be responsible for such a startling effect. 
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1. Introduction 

Experimental measurements of d-Au collisions at the 
Relativistic Heavy Ion Collider (RHIC) and p-Pb colli¬ 
sions at the Large Hadron Collider (LHC) have shown a 
peculiar centrality dependent effect on the nuclear mod¬ 
ification factor {Rpjd-A) for jets and leading hadrons as 
shown in Ref. m There is an observed enhancement 
in peripheral events and a result consistent with unity 
or further suppressed for central events as seen in Refs. 
O, lO, and P4l. This is in conflict with our naive expec¬ 
tation that peripheral events should behave as p-p events 
and that any deviation from unity would be more pro¬ 
nounced in central events due to nuclear effects. 

We present an explanation for these results as well 
as a test of the effects of this explanation on Rpjd-A for 
comparison to these results. In a given p/d - A event 


where there is a jet, the nucleon’s parton distribution 
function (PDF) must be in a state where there is at least 
one hard parton. However, when in this state, there is a 
reduced amount of energy available to produce soft par- 
tons. This depletion of soft partons results in a lower 
number of soft particles being produced for centrality 
binning. In addition, at these ultra-relativistic energies, 
the PDF cannot change (fluctuate) between subsequent 
nucleon-nucleon collisions. This produces a net effect 
that, in an event with a jet, lower numbers of soft parti¬ 
cles are produced than what would have been produced 
if the event did not have a jet. This causes the event to 
be binned as a more peripheral event than it would had 
otherwise been. 

To study this, we setup a Monte-Carlo event genera¬ 
tor for p-Pb and d-Au events. The program starts by 
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generating a heavier nucleus (Pb or Au) and deuteron 
for d - Au events. These are handed off to an impact 
monte carlo that determines preliminary event informa¬ 
tion to be fed into PYTHIA. In order to properly use 
PYTHIA for this analysis, PYTHIA had to be modi¬ 
fied in order to incorporate event-by-event shadowing, 
in addition to a modification to prevent PDF resampling 
during the event. PYTHIA observables were then used 
to determine Ppp^ or R^au respectively to compare to 
experiment. In the remainder of these proceedings, we 
will describe each of these components. We will con¬ 
clude with a discussion of the results from our simula¬ 
tions and an outlook for future work. 

2. Monte-Carlo Event Generation 

The nuclear Monte-Carlo was used to generate nuclei 
based upon an input density distribution. To generate 
heavy nuclei for our analysis {Au and Ph) we used a 
Woods-Saxon distribution 


the colliding nuclei to shift in angle relative to one an¬ 
other. The second (lighter) nucleus was adjusted in posi¬ 
tion by the impact parameter and nucleon-nucleon colli¬ 
sions were counted by determining how many nucleons 
in the heavier nucleus were within a distance of the di¬ 
ameter of a proton from the colliding nucleon. If this 
particular sampling did not result in at least a single 
nucleon-nucleon collision, the Monte-Carlo would re¬ 
set and choose a new impact parameter, repeating the 
previous procedure until a collision was generated. At 
this stage, relative probabilities of proton and neutron 
participation were evaluated, and the nucleon-nucleon 
collision types were determined (eg. p - n, p - p, etc.). 

In order to generate events and event observables, the 
PYTHIA event generator was used i). While it is not 
natively capable of full heavy-ion events, it is possible 
to modify PYTHIA to suit the purposes of this analysis. 
These modifications to this program were made so that 
it could perform event-by-event nuclear shadowing, as 
well as alterations that would allow us to freeze the PDF. 

To implement nuclear shadowing, PYTHIA was al¬ 
tered to incorporate shadowing on an event-by-event ba¬ 
sis. This was accomplished by modifying the PDF calls 
with a factor of 


for Au, the parameters used were a = 0.535 fm, R = 
6.38 fm and for Pb the parameters used were a = 0.546 
fm, R = 6.62 fm) and for the lighter deuteron we used a 
distribution based upon the Hulthen wavefunction 


with parameters a = 0.228 /m“\ b = 1.18 fm~^. The 
parameters and functions used were obtained from Ref. 
Q . To account for the short-distance nucleon repulsion 
force, an exculsion volume for subsequently generated 
nucleons was used; this was set to be the diameter of 
a proton, 1.12 fm. The Woods-Saxon distribution was 
sampled out to 13.3 fm and the Hulthen distribution was 
sampled out to 38.9 fm to ensure that the rare fiuctua- 
tions of nucleons well outside the central dense nucleus 
were accounted for. The large value for the Hulthen dis¬ 
tribution was used to ensure that the long tails of distri¬ 
bution were captured with an error less than 3.67 • 10“^^ 
The generated nuclei were then fed into the impact 
Monte Carlo. The nuclei are projected onto the plane 
perpendicular to their momentum and an impact pa¬ 
rameter was chosen for the event. This was performed 
by sampling an distribution out to 10 fm; distances 
greater than this required more than ten million gener¬ 
ated events to provide even a single collision. In ad¬ 
dition, a random 6 was chosen from 0 to In to allow 


fix) = 1 + (^(v) - 1). (Ncoii/ < Ncoii >). (3) 

where 5'(v) is the nuclear shadowing function from Ref. 
m. This will return the average value of the shadowing 
over a large number of events, but incorporates event- 
by-event fluctuations. The number of collisions was 
used to adjust the shadowing. 

To incorporate the fact that the PDF of an impact¬ 
ing nucleon does not resample during an event, in¬ 
cluding events where a hard parton was produced, the 
event was generated as a single impact of a proton on 
a ’’superproton”. This is in constrast to the more typi¬ 
cal procedure where the event is treated as the sum of 
n nucleon-nucleon collisions. This was accomplished 
by taking the number of nucleons and multiplying the 
PDF calls. By taking into account the number of neu¬ 
trons and protons (determined earlier during the impact 
Monte Carlo stage) and multiplying the quark distribu¬ 
tion calls within the PDF calls, preservation of isospin 
can be maintained. As an example, the up quark distri¬ 
bution within the superproton will look like 

F^Ax) = Np-F,{x) + N„-Fdx). (4) 

for Np protons and neutrons and neglecting shadow¬ 
ing. However, because PYTHIA treats the superproton 
as only having charge -r 1, net charge for the entire event 
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will not be generated correctly. At the energies of inter¬ 
est however, this does not effect the number of charged 
particles generated. 


3. Results 

Before comparing to experiment we constructed a 
plot to quantify the effects of ’’bin-shifting”; where 
events in one centrality class as determined by Ntm were 
shifted to another centrality class if binned by Nchg- 
Generated events were tracked, and a plot of this bin- 
shifting was constructed as a function of the pr of the 
leading from the simulation for RHIC data as shown 
in Figure 
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Figure 2: R^au for Monte Carlo Event Generator for 
minimum bias compared to PHENIX results 
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Figure 1: Percentage change of events in a given cen¬ 
trality bin due to ’’bin-shifting” as a function of pr for 
d - Au events 
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To compare to PHENIX R^au data O, tt^s at y < 
|0.35| and Nchg for each event were recorded. While 
PHENIX determined centrality from Nchg in the forward 
detector El, we determined centrality from Nchg over 
the entire event. This was necessary due to low statis¬ 
tics. However this did not cause the minimum bias RdAu 
from generated events to differ from experiment. For 
each centrality bin, dNjdppdy was calculated. Using a 
pion spectrum from pp events generated by PYTHIA, 
^dAu was calculated. This is shown in Figure and a 
comparison plot for RdAu where centrality was deter¬ 
mined with Nhin is shown in Figure The minimum 
bias RdAu was also plotted to demonstrate the veracity 
of this event generator in the absence of any centrality 
dependent effects; this is shown in Figure 

A similar procedure was performed for p - Pb data 
from ATLAS, noting that there are two sets of data we 
wish to compare. The first being Rcp for inclusive jets 


Figure 3: RdAu for Monte Carlo Event Generator with 
centrality determined by Nchg compared to PHENIX re¬ 
sults 


El and the second was Rcp for charged hadron produc¬ 
tion (H. Experimentally, binning was performed by Ep 
in the pseudorapidity interval 3.1 < p < 4.9 in the di¬ 
rection of the lead beam El for both data sets. Our cen¬ 
trality binning was performed by finding the total en¬ 
ergy deposited outside of the two highest pp jets in the 
event or by collecting all charged particles over the en¬ 
tire event that carried less than 10 GeV of pp. 

For jet Rcp, Fast Jet was used to determine jet prop¬ 
erties for each event. The anti-kT algorithm was used 
with R=0.4. The two highest pp jets were recorded for 
each event, along with the total Ep outside of the top 
two jets. Using these jets, dNjdppdy was calculated for 
the 0-10% and 60-90% centrality bins which were di¬ 
vided to obtain Rcp. The procedure was repeated for the 
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Figure 4: R^au foi* Monte Carlo Event Generator with 
centrality determined by Ntm compared to PHENIX re¬ 
sults 


charged hadron data, except centrality was determined 
by Nchg- These are shown in Figures]^ and 



Figure 5: Rcp from Monte Carlo Event Generator com¬ 
pared to ATLAS results for charged hadrons 


4. Conclusion 

In conclusion, we see that the centrality dependent ef¬ 
fect on RdAu of this event generator when binned by Nchg 
is consistent with PHENIX experimental results. In ad¬ 
dition, when binning by Ntm, we obtain the naive expec¬ 
tation that RdAu for peripheral events is consistent with 
unity and that for central events we see a slight enhance¬ 
ment due to nuclear shadowing. For LHC p-Pb events, 
we find agreement between jet Rcp, however, there is a 



Figure 6\ Rcp from Monte Carlo Event Generator com¬ 
pared to ATLAS results for jets 


discrepancy for charged hadron Rcp. This is likely due 
to the enhanced color connections in the superproton. 

The core of this event generator relies on the fact that 
generation of a hard parton to create a jet results in the 
suppression of the production of soft particles; when a 
nucleon’s PDF is in a state where there is a hard par- 
ton, it depletes the low x portion of the PDF. Since bin¬ 
ning of events is performed by some measure of soft 
particle production, this leads to events with jets be¬ 
ing shifted to more peripheral bins than those events 
without a jet. This was incorporated by ’’freezing” 
the PDF with PYTHIA; treating the event as a proton- 
superproton event invokes PYTHIA s energy conserva¬ 
tion routines over the entire event. With this event gen¬ 
erator, we simulate the observed centrality dependant 
effects on Rpjd-A and Rcp for d - Au and p - Pb events 
as a result of the this ’’frozen PDF” effect. 
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